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ABSTRACT 


Since  the  availability  of  analytic  solutions  of  the  Schrodinger 
equation  are  limited  to  the  most  simple  atoms  and  molecules,  the  necessity 
for  methods  of  approximation  has  been  keenly  felt,  one  such  method  being  the 
use  of  the  linear  variational  method.  However  if  this  method  is  to  be  em¬ 
ployed,  use  has  to  be  made  of  an  initial  trial  function,  the  form  of  which 
is  somewhat  arbitrary,  and  may  depend  upon  any  physical  information 
available  or  chemical  intuition.  Because  of  their  ressemblance  to  the  avail¬ 
able  analytic  solutions,  the  Slater  type  functions  have  found  extensive  appli¬ 
cations;  but  unfortunately  the  integrals  involved  are  usually  rather  difficult 
to  evaluate. 

This  work  describes  the  use  of  Gaussian  functions  with 
which  the  evaluation  of  integrals  is  quite  simplified.  It  was  the  purpose 
of  this  research  to  investigate  the  advantages  and  disadvantages  of  these 
functions  with  respect  to  the  energy  obtained,  the  electronic  distribution 
and  the  derived  results.  For  this  purpose  the  hydrogen  molecular  ion 
was  chosen  as  the  testing  ground  mainly  because  of  the  availability  of  the 
exact  solution  and  a  number  of  theoretical  approximations. 

Three  types  of  Gaussian  function  were  used  in  the  form  of 
linear  combinations,  namely  s  type,  p  type  and  a  form  involving  a  weighting 
on  the  Z  axis;  the  latter  two  were  used  with  an  idea  of  making  the  electronic 
distribution  more  favorable  in  the  internuclear  region.  The  functions  mentioned 
above  have  been  referred  to  in  the  text  as  the  first,  third  and  second  trial 


function.  It  was  found  that  although  the  energies  obtained  frojn  the  first  type 
of  function,  using  a  combination  of  three  Gaussian  functions  (which  is  a 
number  that  can  be  readily  handled  by  the  IBM  1620)  were  quite  good,  the 
electronic  distribution  was  very  unfavorable  chiefly  because  of  the  presence 
of  rather  flat  maxima  at  the  nucleii,  as  against  the  correct  sharp  appearance 
given  by  the  exact  solution.  The  form  improves  with  a  larger  basis*  set; 
but  a  very  large  set  involves  a  very  laborious  computation  and  renders  the 
problem  impracticable  from  the  standpoint  of  time.  In  the  second  type  of 
function  there  is  a  small  decrease  in  the  energy,  whereas  thp  numerical 
integration  involved  is  too  time  consuming  to  permit  a  recommendation 
of  this  function  for  any  practical  use.  In  the  third  function,  tjie  introduction 
of  a  p  type  orbital,  does  in  fact,  cause  an  increase  in  the  electron  density 
in  the  middle  region;  but  the  expected  improvements  are  largely  offset  by 
a  rather  poor  representation  of  the  s -  character  about  the  nucleii.  This 
is  caused  by  the  small  number  of  s  type  Gaussian  functions  used  in  order 
to  keep  the  problem  within  practical  limits. 

It  is  concluded  that  it  is  advantageous  from  the  practical 
point  of  view  to  introduce  as  the  third  basis  function,  one  more  s  type  member 
instead  of  a  p  ."type  function. 

In  spite  of  the  poor  description  of  the  electron  distribution 
obtained  by  the  Gaussians, the  derived  results  such  as  the  spectroscopic 
constants  are  quite  good,  casting  doub.tt  on  the  suggestion  for  their  use  as 
criteria  for  the  ‘goodness*'  of  a  function. 
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CHAPTER  1 


INTRODUCTION 


The  Hydrogen  Molecule  Ion 


The  hydrogen  molecular  ion  has  been  the  subject  of 
a  considerable  amount  of  research,  during  the  last  couple  of  decades, 
mainly  because  it  bears  the  same  relation  to  most  molecular  problems, 
as  does  the  hydrogen  atom  to  most  atomic  problems;  like  the  latter 
it  is  a  single  electronic  species,  free  of  interelectronic  terms,  and 
yet  it  has  a  distinctly  molecular  nature.  It  was  decided  therefore  to 
use  this  ion  as  a  testing  ground  for  the  use  of  the  Gaussian  Functions. 

This  molecule  was  discovered  many  years  ago  by 
J.  J.  Thomson  in  cathode  rays,  produced  by  the  bombardment  of 
ordinary  hydrogen  gas  with  electrons.  Since  then  it  has  received 
a  great  deal  of  theoretical  investigation,  of  which  a  short  review 
would  be  worth  while :  - 

Accurate  solution  of  theSchrodinger  equation: - 

The  system  consists  of 

two  nucleii,  and  one  electron.  From  the  Born  -  Oppenheimer 
approximation,  one  can  neglect  the  motion  of  the  nucleii,  and  write 
down  the  Schrodinger  equation  as 


in  which  and  represent  the  distance  of  the  electron  from 

nucleus  A  and  B  respectively.  If  elliptical  co-ordinates  §>y\  and 
W  defined  by  the  equations 
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(with  *A»  as  internuclear  distance)  be  introduced,  then  the 
partial  differential  equation  becomes  separable  into  three 
total,  differential  equations.  Introducing: 
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the  equation  reduces  to: 
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Where  'Yt\  , ^  and  are  undetermined  parameters.  The  familiar 

equation  possesses  acceptable  solutions  for  m  -s-o,  ±]  . .  . 

The  subsequent  procedure  of  solution  consists  of  finding  the 
relation  which  must  exist  between  and  h-  in  order  that  the 

equation  posses  a  satisfactory  solution,  and  then  using  this  relation, 
for  finding  from  the  equation,  characteristic  values  of  A  . 

This  procedure  was  carried  out  by  Burrau  on  the  normal  state 
of  the  hydrogen  molecular  ion,  by  a  numerical  integration  of  equations 
(3)  and  (4).  More  accurate  studies  have  been  made  on  equation  (4) 
by  Hyelleras  (2)  and  on  equation  (3)  by  Jaffe^3)  and  Chakrava  rty 
From  the  definition  of  \  —  Cos#  ,  it  is  evident 
that  the  only  values  this  variable  can  take  up,  lie  between  -1  and^l, 
which  also  happens  to  be  the  limits  of  the  associated  Legendre 
function  pjj^  (Cos  H  ),  hence  taking  advantage  of  the  complete  nature 
of  these  functions  one  can  write. 
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Substitution  of  this  into  equation  (4)  and  simplification  with  the  aid 
of  the  differential  equation,  satisfied  by  the  associated  Legendre 


function,  one  obtains 


(4) 


I  Q f >t--H—-e(e+i)}pJrn/(n)~o  (*) 

J-pmi  L  J 

Elminating  the  factor  T*  by  successive  use  of  the 
ordinary  recursion  formula  (5),  and  upon  introducing  the  result  so 
obtained  into  equation  (5)  a  simple  series  in  the  functions  P  ^  ^  ) 

with  coefficients  independant  of  results.  Because  of  the 
orthogonality  of  these  functions,  their  coefficients  must  vanish 
independently,  in  order  that  the  series  vanish.  This  gives  the 
condition 

xc  -  f  f  {f+fabi)  ,  (e-fo  )(#+  imj  h 

't  Xe+J  -O  (6) 
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The  equation  (6)  may  now  be  considered,  as  a  set  of 
homogenous  equations  in  the  unknown  quantity  ,  with  different 
values  of  ('  .  In  order  that  this  set  may  possess  a  non-trivial 
solution,  the  determinant  must  vanish.  This  gives  thep  a  deter- 
minantal  equation  involving  /\  and  f*L.  ,  from  which  one  can  determine 
the  relation  between  them. 

For  the  normal  state  of  the  system,  with  m=0  and  £ 


even,  the  determinantal  equation  reduces  to 
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obtaining 


Hyellerass  carried  out  the  solution  of  this  equation, 


^  }>  ^  A-  ^  ~  o •  oo  as  13,  7^-€'"q0 oo  e  a^  ^ 


This  gives  the  functional  relation  between  the  separation  constant  r\ 
and  A  as  derived  from  the  ^  equation,  the  next  step  would  be  to 
do  the  same  with  the  ^  equation  and  then  equate  them.  The  latter 
relation  has  been  derived  both  by  Jaffe  and  by  Chakravarty. 

If  in  equation  (3)  2D  is  replaced  by  B,  and  M'  by  -A, 
the  equation  can  be  readily  rearranged  as 
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setting 


d'rt %  and 
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P 

and  taking  (l-x)'o-  for  A  one  obtains  after  simplification 
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Taking  ^  (which  is  as  yet  undefined)  to  be  equal  to 
the  equation  can  be  further  simplified 

m  -a  fn  )i  x1 + Cfo/v*  pn*  <)  ax— 


-  o- 


(9) 


where 
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The  equation  (9)  is  then  an  ordinary  differential  equation  that  can 


be  solved  by  series  integration,  by  setting 
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From  whence  one  gets  the  recursion  formula 
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From  (10)  then 
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Where 
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Hence  from  (12)  as  a  first  approximation  one  has 
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(13) 


and  to  a  second  approximation,  considering  only  the  first  convengents 
one  has 


(13a 
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Hence  substituting  for  K  its  value  as  given  by  (14)  wherever  it 


(8) 


occurs  in  the  second  order  terms  in  (13a)  ie  in  '’’tn  and 

u,  »  Wi 

\  one  obtains  on  simplification 
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This  equation  can  be  simplified  for  the  ground  state  and 
equated  to  (7)  obtained  for  the  part;  with  the  help  of  these  one 
then  proceeds  to  obtain  the  roots  of  \  from  the  definition  of  which, 
equation  (4a),  and  by  taking  suitable  values  of  one  can  readily 

evaluate  the  E  eigenvalue,  which  agrees  remarkably  well  with 
experimental  values. 

The  eigenfunction  is  given  by 

wheye  N  is  the  normalization  constant. 
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where  the  coefficients  ,  £L  . 
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However,  the  wave  function  obtained  above  is  rather  cumbersome,  and 
difficult  to  handle,  and  it  would  be  desi  reable  to  have  an  approximate 
treatment  of  the  molecule. 

Morse  and  Stueckelberg 's  degenerate  perturbation  treatment:- 

One  of  the  earlier  attempts  to  study  the  hydrogen 
molecular  ion  by  an  approximation  method  was  made  by  Morse  and 
Stueckelberg  in  1929^)  using  the  perturbation  method.  These 
workers  took  advantage  of  the  fact,  that  at  least  under  two  conditions, 
the  eigenvalues  and  eigenfunctions  of  this  species  are  known  precisely 
(l)  when  the  internuclear  distance  is  zero,  one  gets  the  He  ion,  which 
has  exactly  the  same  solutions  as  an  hydrogen  atom  carrying  a 
nuclear  charge  two,  and  (2)  when  the  internuclear  distance  is 
infinity,  one  gets  essentially  an  hydrogen  atom  and  a  proton  at  an 
infinite  distance  apart,  the  eigenfunctions  and  eigenvalues  of  the 
of  the  hydrogen  atom  are  of  course  well  known.  In  between  these 
two  extremes,  there  lies  the  hydrogen  molecular  ion  being 
investigated. 

The  change  over  from  the  He  ion  at  a  zero  nuclear 
separation  to  an  hydrogen  atom  and  a  proton  at  infinite  nuclear 
separtion  via  the  H^  lying  in  between,  was  considered  by  Morse  and 
Stueckelberg  in  the  light  of  the  best  co-ordinate  system,  that  can  be 
used  to  describe  a  particular  species  encountered,  eg.  an  He  ion  is 
best  described  in  a  spherical  co-ordinate  system,  while  an  hydrogen 
atom  with  a  proton  at  an  infinite  distance  from  it,  would  on  the  other 
hand  be  best  described  in  parabolic  co-ordinates,  and  in  between  the 


hydrogen  molecular  ion  is  best  suited  for  an  elliptical  system. 

Since  at  a  finite  distance  from  the  nucleus,  each 
quantum  number  corresponds  to  a  node,  use  can  be  made  of  the 
latter  to  keep  track  of  the  former  during  this  process  of 
changing  internuclear  distances. 

Starting  then  from  the  He  ion,  which  might  also 
be  called  the  united  atom,  one  has  the  spherical  axes,  and 

,  with  the  corresponding  quantum  numbers  ^  and  . 

In  the  modern  system  of  nomenclature 

7/1=  7l*j? 

T\  s  "h>i,  -f  T|  ©  "f  “Y\(p  -f  | J 

or  in  other  words  the  total  quantum  number  is  equal  to  the  total 
number  of  nodes  plus  one.  The  energy  is  given  by 

EW  -  V-™-8*  _  R  =  7-^ - -jt 

where  R  is  the  energy  of  the  hydrogen  atom.  The  yvave  function 
here  is  the  hyjdrogenic  wave  function  with^—  2.  Now  suppose 
the  internuclear  distance  is  increased  from  zero  to  a  finite  value 
of  4^  g,  ,  producing  thereby  an  ion.  One  has  in  effect 

moved  from  a  spherical  system  to  an  elliptical  system,  and  use 
can  be  made  of  the  corresponding  co-ordinate  system  defined 


earlier.  By  taking  limits  of  these  co-ordinates  one  finds: 
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It  is  possible  then  to  write 

^6>  "* 

^cp-*  flip 

and  the  total  quantum  number  is  still 


But  of  course  there  is  no  enegy  expression  available. 

When  the  nucleii  are  moved  still  further  apart  to 
infinity,  one  ends  up  with  a  hydrogen  atom  at  A,  proton  at  B,  or  a 
proton  at  A  and  a  hydrogen  atom  at  B.  The  wave  function  describing 
this  system  will  be  a  linear  combination  of  these  two  forms,  i.e. 

where  both  plus  and  minus  signs  have  to  be  taken,  since  the  probability 
is  given  by  the  integral  of  and  the  pres  ence  of  either  of  these 

signs  makes  no  diffferpnce.  If  one  defines  a  set  of  parabolic 
co-ordinates  at  each  center  as 

/to  (ti  tedn)  t  f 

and  the  same  for  the  B  center,  Then  again  by  taking  limits  one 


finds 
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Changes  in  the  various  quantum  numbers  during  this 
process  may  again  be  found  by  an  examination  of  the  nodes.  The 
Yy  remain  unthanged,  all  theT|iget  converted  into  but  since  the 
TV  are  a  family  of  hyperboloids  of  revolution,  it  would  be  expected 

v 

therefore,  when  separation  is  effected,  that  if  "^r^happens  to  be  even, 
one  half  of  these  would  go  with  one  nucleus,  and  the  other  half  with 
the  other  nucleus,  i.e.  r'  q  ?  V  ;  however,  if  Tjyj  is  odd 

then  one  of  the  hyperboloids,  stays  in  the  middle,  and  a  node  appears 
there;  either  side  of  which  the  wave  function  takes  on  positive  and 
negative  signs,  thus  corresponding  to  the  anti- sym me tark  state. 

One  has  now  a  system  of  a  proton  and  an  hydrogen  atom,  with 
total  mumber  of  nodes  equal  to  (  ^^-{-7^  )  or 

(  A ^  4-Ty -fZ.  )•  The  energy  being  that  of  the  hydrogen  atorn^ 
the  expression  for  which  is  well  known,  and  so  are  the  corresponding 
eigenfunctions , 

<  .  * 

E  i-  — — — —  7 - - - 

t^iyn  xn^  ij2*  >  -f-’lf + Xtof&r 


One  therefore  has  a  knowledge  of  the  system,  at  the 
two  extreme  ends,  and  the  corresponding  wave  functions  exactly. 
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The  next  step  is  to  move  a  bit  from  either  extremes,  immediately 
potential  energy  terms  of  the  type  V  T  — ^7^  begin 

to  appear  in  the  Hamiltonian,  and  one  then  usesthe  extreme  wave 
functions  for  a  degenerate  perturbation  calculation,  aid  thus  obtains 
a  few  more  points  close  to  the  ends  (one  cannot  carry  out  these 
calculations  too  far  in  because  the  perturbation  terms  become 
very  large).  Then  with  the  extreme  ends,  and  the  few  other  points 
calculated,  one  can  draw  a  smooth  curve,  which  can  be  used  to 
calcute  the  energy  at  any  internuclear  distance. 

Using  this  method  Morse  and  Stueckelberg  found  a  very 
close  similarity  with  the  state  obtained  by  Burrau.  They  also 

found  that  minima  appear  only  in  the  ,  3d  &  and  4f  6  states. 

Therefore  to  this  degree  of  approximation,  the  molecule  H^*  is  only 


stable  in  these  states. 


Gullimi  n  and  Zeners  Variational  Study:- 

About  the  same  time  as 

Mose  and  Stueckelber'g  carried  out  the  above  analysis,  Gullimin  and 
Zener  (7),  carried  out  a  variational  study  of  the  molecule,  and  used 

the  same  line  of  argu  ments  to  arrive  at  an  intial  guess  function. 

> 

These  argu-  ments  are  as  follows: 

(l)  The  differential  equation  describing  the  hydrogen  molecular  ion, 
is  separable  into  elliptic  co-ordinates  ^  and  V  suggesting 
thereby  the  possibility  of  taking  for  the  eigenfunction  a  product  of 


two  co-Q\rdinate|functions . 


. 


(14) 


^  LC§)M^n) 


The  absence  of  degeneracy  in  the  ground  state  function  strengthens 
the  assumption,  as  it  precludes  the  neccessity  of  having  to  take 
linear  combination  of  all  degenerate  states,  and  thereby  making 
the  above  factorization  impossible 

(2)  The  equation  being  even  in  ^  ,  undergoes  no  transformation 
when  operated  upon  by  the  inversion  operator  or  f  (n)=t(  -\) 
enabling  one  to  take  either  the  positive  or  the  negative  function 
with  equally  good  assurance  of  describing  the  species,  but  the  com¬ 
plete  absence  of  nodes  in  the  ground  state,  favors  a  positive  function. 

(3)  Reduction  of  the  internuclear  distance  to  zero  results  in  the 

■f 

production  of  a  He  ion,  having  a  hydrogenic  wave  furction,  with  a 
nuclear  char  ge  of  2 


where  r  is  the  distance  of  the  electron  from  the  middle  of  the 
nuclear  axis.  For  an  infinite  separation  however,  one  gets  a  hydrogen 
atom  at  A  and  a  proton  at  B  or  vice  versa,  which  obliges  one  to  take 
a  linear  combination  of  the  type 


1  (  e'^  +  e  r6  ) 

4\fn 
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Conversion  to  elliptical  co-ordinates  defined  earlier,  gives 

c (oust)  e'/J ^ * ( e "4 e 


Here  one  can  say  that 


Taking  limits  of  this  function,  one  finds 


-L 

**”>  \a 


=  %Hr 

-  o 


AM 

Pn 
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at  one  center 


/»<& 


- 


at  second  center 


••  *«cCe%e*) 


for 


Hence  the  wave  function  degenerates  into  one  describing  a 
-L 

He1  ion  and  one  describing  a  hydrogen  atom  and  a  proton,  for 
zero  and  infinite  nuclear  separations  respectively,  therefore 
satisfying  the  physical  conditionsrequired.  This  guessed  function 
is  then  substituted  into  the  expression  f  ^  and  the  energy 

minimized  with  respect  to  the#"  and  B  parameters  (results  are 
tabulated  later) 

More  recent  perturbation  calculations:- 

Although  no  new  numerical 
information  regarding  the  hydrogen  molecular  ion  can  be  obtained 


'  ' 


* 
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by  treating  the  wave  equation  by  perturbation  methods,  never  the, 
lessimportant  information  regarding  the  physical  form  of  the  ion, 
was  obtained  by  using  such  methods,  and  they  also  lead  to  better 
variational  guess  functions.  The  perturbation  methods  can  be 
applied  to  many  systems  for  which  the  wave  equation  cannot  be 
solved  directly,  and  study  of  these  methods  using  the  hydrogen 
molecular  ion  as  a  test  species  was  considered  important  from  a 
comparative  point  of  view.  Treatment  along  these  lines  was  first 

carried  out  by  Pauling  0)(^) 

One  proceeds  then  by  making  reasonable  giiesses  at 
thewave  function  for  H^h  calculating  the  dependence  of  the  energy 
on  the  internuclear  separation  in  each  case.  If  the  energy 

passes  through  a  minim<//j»  at  some  value  of  ,  the  formation 

of  a  stable  molecule  is  indicated. 

Assuming  as  a  first  attempt  that  the  wave  function 
is  simply  that  of  an  electron  localized  in  a  IS  -orbital  about  one  of 
the  nucleii  say  A,  one  has 

.  i  -  111.  - 

tf  =  wis*  -  tt  e 

The  energy  corresponding  to  this  function  is 


E  _  j  H  /  / 

<  uts*  I  K/% y 


(17) 


since  is  normalized,  the  denominator  cbn  be  set  to  unity 

E  =.  [-  ii—  V  X-  /  U,*,  )  +<>*  [ \B  [u«*> 

'Cls«  I  /U' Sfl) 

The  first  i  ntegral  here  is  the  energy  of  the  hydrogen  atom 
R~  -l/Z  atomic  units,  the  second  integral  is  I/^ab  the 
application  of  the  Borm -Oppenheimer  approximation,  while  the 
third  integral  can  be  evaluated  by  using  an  electrostatic  model 
ie  of  having  a  unit  cha^g.  at  a  distance  ^^from  the  center  of  an 
electron  cloud  whose  charge  density  is  I  ^/S 

Summing  all  these  integrals  together  one  obtains 


It  is  evident  that  a  reduction  of  „  to  zero  would 
result  in  an  infinite  energy  value,  which  in  physical  terms  could 
be  considered  as  being  caused  by  the  penetration  of  a  proton 
into  a  negatively  changed  atmosphere,  carrying  another  proton 
in,  its  middle,  and  thereby  displaying  a  completelack  of  stability 

-4. 

of  the  H^' molecule.  This  wave  function  then  does  not  describe  the 
situation  accurately  at  all. 

Two  of  the  essential  features  would  be  (l)  the  possibility 
that  the  electron  can  pass  back  and  forth  from  one  nucleus  to  the 


(18) 


other  and  (2)  the  possibility  that  the  electron  will  tend  to  spend 
much  of  its  time  in  the  region  between  the  two  nucleii,  where  its 
potential  is  low.  Hence,  forming  a  linear  combination  of  j 

and  one  obtains  the  secular  equation 


0 


where 

haa  ~  <CuiSa  j  hj 

HAB  -  \  /ti)  V‘S$&y 

A  Ui$*  f  M 

&  represents  lack  of  orthogonality  of  and  ^  Since 

the  two  nucleii  are  identical  one  has  HAA  ^BB  anc^  since  H  is 
hermitiam  and  real  HAB==HBA,  one  can  obtain  two  solutions  of  the 
secular  determinant 


E„  H’AA  -t"  HAB 

S=  - l-f 

E  A  HAA  '  HAB 

1  -  A 

These  correspond  to  the  wave  functions 

1  ( ) 

\Z  -r.Z&T 

1  (u^-a^) 

-  2  a 
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where  S  and  A  represent  the  words  symmetric  and  antisymmetric 
wave  functions  respectively.  The  wave  function^  is  symmetric  in  the 
positional  co-ordinates  of  the  two  nucleii  A  and  B.  This  function  places 
the  electron  with  equal  probability  in  the  vicinity  of  each  of  the  two 
nucleii.  Thus  in  the  plane  midway  between  the  two  protons  one  has 

,  and  if  the  two  nucleii  are  far  apart,  /^| 

The  electron  density  in  the  midplane  of  the  molecule  therefore 
has  approximately  double  the  value  that  would  be  produced  in  the 
same  region  by  a  simple/S  orbital.  If  one  parries  out  the  integrations 
involved  in  the  expression  for  >  one  obtains 


5,  ^ 


j_  ,  (^b  ± 0  wv'i  £  -ai 


is  very  large,  it  gives  an  energy  of  -l/2  which  is  the 
energy  of  an  isolated  hydrogen  atom.  As  /*^gis  decresed,  the 
energy  passes  through  a  minirrun.  The  calculated  equilibrium 
interfiuclear  distance^  and  the  dissociation  energy  are  comparable  in 
order  of  magnitude  with  the  experimentally  observed  values. 

By  allowing  the  electron  to  move  in  the  vicinity  of 
both  nucleii  one  has  in  effect  introduced  the  resonance  integral 

Ha„,  which  contributes  to  the  formation  of  a  stable  bond. 

AB 

The  chief  stabilizing  factor  in  this  wave  function  is  a 
reduction  of  the  potential  energy  when  the  electron  is  in  between  the 


two  nucleii. 


. 
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binding  energy  was  quite  a  bit  lower  than  would  be  expected  from 
experimental  studies.  Recently  however  Walsh  and  Matsen  (9) 
have  suggested,  that  owing  to  the  presence  of  strong  charges  at 
the  points  A  and  B,  there  will  be  a  strong  force  tending  to  make  the 
electron  prefer  the  direction  of  the  internuclear  axis,  and  they 
describe  the  low  bindings ene rgy  obtained  above  to  the  fact  that  the 
wave  function  lacks  a  polarization  in  this  direction,  and  hence  they 
considered  hundred  percent  polarization,  by  assuming  the  inter¬ 
nuclear  axes  to  lie  in  X-direction,  and  carrying  out  Pauling's 
treatment  in  this  set  up,  using  one  dimensional  wave  furctions.  One 
can  calculate  the  exchange  integral;but  not  the  Coulomb  integral. 
Using  the  exchangeintegral  so  obtained  and  the  three  dimensional 
Coulomb  integral,  these  workers  report  a  37%  higher  binding  energy 
than  the  exact  value. 

Variational  Studies  of  the  H-J~  Molecule:  - 

The  type  of  wave  function 

that  has  been  used  in  the  Pauling  treatment,  does  not  give  one 
accurate  description  of  the  species.  This  discrepancy  is  not 
surprising  in  view  of  the  crudeness  of  the  function.  The  electronic 
cloud  as  given  by  the  wave  function,  is  of  the  same  shape  as  that 
obtained  by  the  accurate  solution,  but  it  is  too  small  in  the  region 
between  the  two  protons,  and  toolarge  elsewhere. 


The  simplest  improvement  is  the  introduction  of  a 


' 


scale  factor?  into  the  hydrogenic  15  functions 


this  give£3a  parameter  ^  ,  with  respect  to  which  minimization  can 
be  effected.  This  wave  function  was  used  by  Finkelstein  and 
Horowitz  (10),  with  very  marked  improvements,  as  will  be  seen 
from  the  table  of  comparison  to  follow  (end  of  Chapter  IV). 

4- 

The  value  of^  that  corresponds  to  the  stable 
molecule  is  1.228,  in  the  Pauling  function  ^jris  given  the  value 
unity.  The  role  of  this  ^  factor  is  to  bind  the  electron  cloud  up 
more  strongly,  in  the  region  between  the  nucleii.  At  other  inter- 
nuclear  distances  the  optimum  value  of  2.  has  different  values.  If 
the  internuclear  distance  is  zero,  one  has  an  helium  ion  in  its 
ground  state,  for  which  the  correct  wave  function  is  a  /S  orbite.1 
with  a  nuclear  charge  of  2,  ie  ^  2.  On  the  other  hand,  if  the  two 

nucleii  are  infinitely  for  apart  the  wave  function  must  reduce  to 
that  of  a  normal  hydrogen  atom  with  1.  The  scale  factor  may 
be  regarded  as  an  effective  nuclear  charge  whose  value  increases, 
as  internuclear  distance  is  decreased,  because  of  the  electrical 
attractions  that  the  nucleii  exert  on  the  electron. 

Further  improvement  of  the  function  was  effected 
by  replacing  the  pure  IS  atomic  orbitals,  by  hybridized  atomic 
orbitals  of  the  type  (  )  where  C  is  a  variable  parameter 


. 
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(to  be  determined  by  minimizing  the  energy)  and  ^g^is  a  hydrogen 


like  2 p  atomic  orbital,  with  axis  lying  along  the  molecular  axis. 

The  use  of  th^s  function  lies  in  the  fact  that  the  2j^  orbital  has  a 
positive  lobe,  which  can  be  used  to  increase  the  electron  density 
in  the  internuclear  region,  where  the  electron  wants  to  be,  at  the 
sametime  the  presence  of  a  node  at  the  ends  of  the  bond 
reduces  the  probaba.bility  of  finding  the  electron  there’,  this 
happens  to  be  the  place  where  the  electron  wants  to  be  least.  The 
2p^  function  can  be  said  to  make  possible  the  mutual  polarization  of 
the  two  atomic  orbitals. 


(H) 


Dickinson 


has  calculated  the  ene  rgy  corresponding 


to  the  hybridized  function. 


This  function  involves  three  parameters  "i.  ,  C,  and  ,  where  ^ 


and  are  the  angles  between  the  internuclear  axis,  and  the 

lines  joining  the  nucleii  to  the  electron.  This  wave  function  does 
indeed  as  expected  move  the  electron  into  the  middle  and  away 
from  the  two  ends. 


increase  the  electronic  cloud  in  the  middle;  He  used  the  same 
kind  of  wave  function  as  used  by  Pauling,  except  he  centered 
them  at  a  distance^,  inside  from  each  nucleus. 


(24) 


where  and  are  the  distances  of  the  electron  as  measured  from  the 

displaced  centers.  Using  this  wave  function  Hurley  carried  out  a  minimization 

\ 

with  respect  to  X  and 't  ,  and  found  the  electronic  cloud  almost  identical 
with  the  exact  theoretical  calculation.  Unfortunately  this  wave  function 
is  not  very  accurate  close  to  the  nucleii,  at  the  same  time  some  of  thel 
integrals  are  very  difficult  to  evaluate  owing  to  the  presence  of  the  factor 

X 

(13) 

James  has  used  a  very  simple  function  of  the  type 

{¥  -  Jfe*  L-*  ^  ^ f  Wj 

The  quantity  zip  has  its  minimum  value  along  the  straight  line  segment 
drawn  between  the  two  nucleii,  and  it  is  constant  at  all  points  along  this  line. 
Therefore  the  above  function  is  constant,  and  has  its  maximum  value  along 
the  axis  of  the  bond;  obviously  it  over- does  the  concentration  of  electrons 
in  the  bond.  Furthermore  it  does  not  lead  to  the  correct  wave  function  in 
the  limit  of  large  internuclear  distance.  James  found  that  this  simple  function 
is  as  good  as  the  Finkelstein  -  Horowitz  LCAO  Junction  if  the  nucleii  are 
not  too  far  apart. 

The  United  Atom  Treatment:- 

Since  the  LCAO  approximation  when  applied  to  H2  molecule 

is  simpler  than  the  calculations  of  Burrau  and  of  Charkravarty  but  the 

(19) 

results  are  not  as  accurate,  Matsen  carried  out  a  simple  treatment 
(called  the  united  atom  treatment)  which  gives  results  in  good  agreement 
with  the  accurate  calculations.  This  method  consists  in  representing  the 
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states  of  by  atomic  orbitals  centered  at  the  middle  of  the  bond. 

Matsen  did  a  perturbation  calculation  with  silch  atomic 
orbitals.  The  Hamiltonian  is  written  as 


where  is  the  distance  between  the  electron  and  the  midpoint  between 
nucleii  A  and  B.  The  first  two  terms  in  the  Hamiltonian  make  up  H  and 
the  last  three  H/  .  The  zero  order  functions  are  the,  hydrogen  like  wave 
functions  with  charged.  .  For  these  states  to  which  first  order  non-degenerate 
perturbation  theory  can  be  applied,  the  energy  i  s  found  with  the  help  of  the 
virial  theorem 

I  I  h**/ 

where  1  -  2,  this  produces  very  good  results  at  short  internuclear  distances. 

Mats©n  found  that  by  carrying  out  a  variational  treatment  with  respect 

to  the  t  parameter,  one  could  improve  the  energy. 

The  Use  of  Gaussian  Functions 

On  examination  of  the  wave  functions  that  have  been  used 

in  all  the  treatments  given  above,  one  finds  that  essentially  all  of  these 

-Ch. 

are  of  the  form  Ne  or  linear  combintations  thereof,  these  generally  result 

(15) 

in  integrals  which  are  very  difficult  to  evaluate.  Boys  for  the  first  time 
demonstrated  a  method  of  successive  approximation  by  which  the  stationary 
state  electronic  wave  function  for  any  configuration  of  atoms  can  be  calculated 
to  any  desired  degree  of  accuracy  by  inclusion  of  a  sufficient  number  of 
terms.  This  method  does  not  depend  on  any  numerical  integration. 

Boys  for  the  first  time  showed  that  all  integrals  between  Gaussian 


' 
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functions  and  Schiodinger  operators  can  be  evaluated  very  readily, 
the  most  complicated  being  that  of  the  electronic  interactions 
between  one  product  of  two  Gaussians  on  different  centers,  with 
another  product  of  two  other  Gaussians.  These  integrals  and 
the  simpler  ones  required  were  all  evaluated  explicitly  in  this 
paper. 

Boys- proposed  that  these  Gaussian  functions  in  the 

most  general  case  be  compounded  into  Slater  determinants, 

followed  by  evaluation  of  integrals  and  minimization  of  the  energy 

by  the  Ritz  variational  method,  which  consists  in  evaluating  the 

integrals  with  a  trial  function  (  in  this  case  the  Gaussians)  that 

depend  on  a  number  of  parameters,  and  varying  these  paraimeters 

until  the  expectation  value  of  the  energy  is  a  tninimum.  The  result 

is  an  upper  limit  to  the  ground  state  energy  of  the  system,  which 

is  likely  to  be  close  if  the  form  of  the  trial  function  ressembles 

(16) 

that  of  the  eigenfunction 

Several  approximate  and  non-converging  methods 
of  evaluating  wave  functions  of  molecules  have  been  widely  used  for 
making  deductions  on  the  general  nature  of  molecules.  The  best 
known  of  these  are  the  molecular  orbital  and  the  localized  bond 
methods.  It  has,  however,  never  been  possible  to  make  a  prior 
predictions  by  such  procedures  in  cases  other  than  the  very 
simplest,  since  the  requi  site  integrals  could  not  be  evaluated,  and  had 
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to  be  estimated  empirically.  However,  if  the  orbitals  used  in 
any  of  these  cases  were  taken  to  be  linear  combinations  of 
functions  derived  from  Gaussian  functions  all  the  required  integrals 
could  be  evaluated  explicitly.  The  analysis  given  by  Boys 
therefore  provides  for  the  first  time  a  method  of  applying 
quantitatively,  and  in  a  fundamental  way,  the  approximate 
methods  of  the  molecular  orbital  and  localized  bond  types. 

Prior  to  the  announcement  of  Boys,  the  only  cases 
where  explicit  methods  were  available  without  numerical  integration, 
were  for  the  two  electron  problem  as  He  and  H^  but  no  converging 
methods  were  available  for  more  general  molecules. 

One  can  define  the  complete  class  of  Gaussian 
functions  to  consist  of  all  functions,  and  of  all  linear  combinations 

fi  y| 

of  functions  of  the  type  X  $  ^  6" a  )  where  1,  m,  n,  have 

any  integral  values,  and  "a"  has  any  positive  value,  and  one 
has  the  general  notation. 


Ka  ~  X“  A-*  pXc- 


& 


--  +(/* 4  2  * 


and  where  A  ,  A  ,  A  have  any  real  values.  By  taking  these 

*  &  - 

functions  and  multiplying  them  byo(  and  i.e.  spin  up  or  spin  down 
functions  one  can  very  readily  construct  the  Slater  determinants,  and 
the  Schrodinger  integrals  between  them  can  be  readily  evaluated. 


the  only  step  of  difficulty  being  the  evaluation  of  the  integral  of 
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the  type 


<N<?»>  <W*/&>  Q*IvlP> 

where  g^,  g^,  g^  and  g^  are  members  of  the  complete  class  of 

Gaussian  functions.  K,  V,  and  M  correspond  to  kinetic  and 
potential  energy  operators  respectively.  If  the  integrals  belonging 

to  this  class  of  Gaussians  are  evaluated  using  simple  functions  of 

_  3- 

the  type  C  “  then  one  can  easily  derive  the  integrals  for  the  general 
class  by  differentiation  with  respect  to  the  parameters  in  the 
integral  formulae,  followed  by  taking  appropriate  linear 


combinations. 


One  might  then  say  that  the  introduction  of  these 


Gaussian  functions,  has  for  the  first  time  provided  one  with  a  quan¬ 
titative  method  of  studying  any  nuclear  frame  work,  and  there  is 
no  need  to  resort  to  very  complicated  numerical  integrations,  as 
one  had  to  in  the  case  of  the  more  conventional  wave  functions 

i 

that  we  re  used  in  the  past. 

Following  Boys  announcement  of  the  Gaussian  integrals 

(17) 

these  wave  functions  were  used  very  effectively  by  Meckler 


for  calculating  the  oxygen  electronic  levels. 

In  a  later  communication  Boys 


(18) 


showed  that 


the  Gaussians  can  be  suitably  modified  to  include  correlation 
terms  of  the  type  —  distance  between  the  ^  and 


.  . .  no 


9V 


(29) 


electrons)  and  still  the  integrals  reduce  to  a  product  of 

elementary  functions  and  error  functions. 

(19) 

Finally  Singer  has  shown  by  a  different 
method  that  the  potential  energy  integrals  associated  with  the 
most  general  kind  of  Gaussian  functions. 


t  =2- Pfe.  **{>(-«> s' 

R  & 

can  be  reduced  to  simple  numerical  quadratures  (  is  any 
positive  definite  quadratic  form  in  the  space  co-ordinates  of  n 
electrons,  is  a  polynomial  in  3m  (  m<^  n)  co-ordinates). 

The  absence  of  any  restriction  on  the  nature  ofQ^, 
leads  to  a  very  interesting  feature  ie  one  can  study  furctions 
of  non-spherical  nature,  since  actually  in  Boys  studies  only 
spherical  Gaussian  could  be  considered.  By  using  Singers 
modifications  it  is  possible  to  consider  Gaussian  with  an 
ellipsoidal  symmetry,  and  thus  introduce  into  the  wave  functions 
a  certain  degree  of  polarization. 

The  present'  investigation  then  was  carried  out  on 
the  hydrogen  molecular  ion  to  study  tpw  spherical,  ellipsoideal 


and  hybridized  Gaussians  could  be  used,  with  a  view  to  obtaining 
their  relative  merits  as  variational  wave  functions 


. 


r 


CHAPTER  II 


GENERAL 


METHODS 


(30) 


Methods  of  Integration:- 

In  the  present  investigation  the  following  functions, 
centered  on  the  individual  nucleii  were  used:- 


(1) 

e 

-  Cx*+f)  - 

(2) 

e 

(3) 

- 

st  C04&  e 

The  method  consisted  of  determining  integrals  of 
the  following  type:- 

The  integrals  for  the  first  type  of  wave  function  were  easily 
obtained  by  straight  forward  integration  between  appropriate  limits, 
as  demonstrated  by  Boys^),  and  resulted  in  products  of  constants 
with  each  other,  with  the  error  function,  and  the  exponential 
function.  The  error  function  was  determined  by  an  extension  of  the 

( 19) 

method  given  by  Hastings  ,  while  the  exponential  function  was 
evaluated  by  multiplication  of  appropriate  tabulated  values. 

The  integrals  for  the  third  type  were  obtained 
quite  readily  from  those  of  the  first  kind,  by  a  simple  process  of 


' 


(31) 


differentiation  as  follows 


Boys 


If  one  defines  the  following,  in  the  manner  done  by 


+  2#  with  =  iT.-flu)*' etc. 

n&  -  Xp  with  CX-Bjcf  etc. 


then:- 


In  the  case  of  the  second  type  of  wave  function,  the 
overlap  and  the  kinetic  energy  integrals  can  be  obtained  by  straight 
forward  integration;  however  the  potential  energy  integrals,  have 
to  be  evaluated,  by  a  slightly  simplified  form  of  the  method  given 
by  Singer 

The  main  value  of  Singer's  technique  lies  in  the  fact 
that  it  can  be  used  for  the  most  generalized  form  of  gaussian 
functions  eg.  any  function  of  the  type 

E  f*  *K“  ®0 

fc  ft 

where  is  any  positive  definite  quadratic  form  in  the  space 


(32) 


co-ordinates  of  electrons,  is  a  polynonial  in  3^yi  co-ordinates. 
Suppose  then  one  is  interested  in  the  evaluation  of  the  potential 
energy  integral 


where  C  is  the  position  of  the  nucleus,  and  A  and  B  are  the  centers  of 
of  two  gaussians  respectively,  one  has  then 


For  a  single  electron  problem  such  as  the  present  one,  let  the 
electronic  and  nuclear  co-ordinates  be  ,  2:  ,  and  Cj(  ,  C ^  ,<ft 

respectively.  Then  from  vectorial  considerations  one  has:- 


where  and^S,  are  the  diagonal  matrices 


d  - 


£  - 


boo 

O  d  O 


) 


refers  to  transpose) 


By  carrying  out  the  matrix-vector  multiplication  in  (l)  one  has 

©  +  Tfltfl 

Fet  o(  p  ~  €,  o.yvd  —  ft  ^  HR 


«v 


♦  > 


<1.  =  Jx-S* 


_  af**  -+e, 


& 


■ 


(33) 


<£  00  <« 


-  *  (.<*.+$_)  x  ceR  £a+<% 

£  dr 


By  defining  the  following: 

o<  -v£_  =  -2 


e«  +  e» 


-  e 


A. 


a,  a|V  +  aB  -  ■x.ZTc-ie.  *  +« 


p  =  d 1  e 


and 


G  -  e  If  £ 


One  has 


-i 


--  fe  +<£) 


-i 


A 


(«f>*+bB^)1'  Cc^+  Jg*/",  C-f  ft* 

»  O  D"1  P  —  _ I-  - - - -  T.  -L  - - - - — 

*'  (o.+  b)  +  CC+ol) 


But  from  the  definition  of  P  one  has 

€*  *  2.  9  . 

^  is/ 

e  =  p  £  c 

£7  =  G  .  p  P  eTd*^  =•  G  -  P  D  P 

o-r  G>  -  P 


(since  D  is  diagonal) 


If  the  vector  between  the  origin  of  the  co-ordinates  and  the 

electron  be  represented  by  JL,  and  that  between  the  origin  and 

— ? 

the  nucleus  by  C,  then  the  vector  between  the  electron  and  the 
nucleus  will  be 


a>c  -  n .  -  c 

-  |  K  -  c  I 


using  the  identity 


cO 


-oo 


one  has 


t  \ 


n ]<u 

60 

0-<*  )z4  Q-  <g  Ci  -  )*■ 


Defining: 

b  a  C  -  P 

One  has 

t  tjtl  P*. 

c>  =  b3+<Y 

c4  =.  bj  ■+  pt 


1*  .  -c| 


,/  T.  „ S- Rr- r*)*4b* - it* l*-P* )  +<**»> +fcL7  , 

rhJ *n  ‘‘ 1^-,’  J' 


. 


(35) 


Substitution  of  this  into  the  operator  part  of  the  integral  under 


investigation  yields, 


-oo 


f  L  oj-fr  £4  dtu-8- 

[(<x4b  +  u2  )C^4-cL  4-u.^  (£43.  -v 


oLc 


Changing  variables,  by  making  the  substitution 


iv-  »-*)*■ 


one  gets, 


ob  ($*— Rk)2* 
o-4  t> 


— r*d  *  +  f  ^ —  jx 


H 


(fl.4b  )b». 


[fo,+b>  +  tfr +*R{CfnMctc^}^J^ 


Under  the  conditions  dictated  by  the  second  form 
of  function,  used  in  this  work,  one  has  (a4*b)  -=(c4d)  -if  (f-fg)  which 
can  be  achieved  by  setting  a=.c  and  bad 

.  «•«**  +£>-f+i  |  »/\J  €  J*»>= 


The  last  integral  in  this  expression  cannot  be  analytically 
evaluted  and  one  is  obliged  to  carry  out  a  numerical  integration. 


■ 


(36) 


During  the  course  of  the  present  investigation  the  Gaussian 
quadrature  method  as  described  by  Hartree  ^^was  used. 

This  method  consists  essentially  of  using  a  finite 
nurrber  of  values  of  the  integrand,  and  may  be  regarded  as  giving  a 
weighted  mean  of  these  values. 


b 

pcx)  Ax  =  U-*)  [ 

a. 

where  the  W  s  are  the  weights  assigned  to  the  values  of  the 

.K 

integrand  at  points/^  . 

Introduction  of  the  new  variable 


J  =  Qlx-O*  b) 


(2) 


(3) 


transforms  the  limits  of  integration  to  -1,  Vl,  and  the  integral 
becomes : 

*“  “I  i 

The  method  next  consists  of  replacing  the  function  <(>  (5)  bV  a 
function  that  can  be  expanded  into  a  series  of  the  type  (2) 

+  1 

f  =  2  (5) 

“I 

Suppose  that  the  value  of  S  is  known  at  (n+l)  points,  then  it  will  be 


(37) 


shown  that  one  can  replace  4>  (t»)  by  a  polynomial  ^  ) 

where  Z  ri+1  is  the  order  of  this  polynomial.  Let  F*n(§  )  be 
a  Legendre  polynomial  of  order  n;  then  from  a  well  known 
orthogonality  relation  one  has 


+1/ 


x*\  w 


(6) 


Dividing  C  §  )  by  gives 

-H?>  *  1.(5)  +  -MIL 


(7) 


where  (^)  and  ^-y\  (^)  are  the  quotient  and  remainder 

polynomials  respectively,  of  orders  n,  which  can  be  proved 


by  induction. 


••  £ 


In 4| 


”  P-yi  +  l  (~tl  )  ^yvC  £,)  +  %■*(  5  ) 


(8) 


Since  the  Legendre  functions  form  a  complete  set,  the  function 
can  be  represented  as  a  linear  combination  of  these 


'n 


(9) 


since  m  is  never  allowed  to  equal  (.n+1)  it  is  obvious  that  the 


integral 


+\ 


f  =o 


(10) 


. 

5 


(38) 


(11) 


Hence  integration  of  both  sides  of  (8)  gives 

+/  +( 

If  one  therefore  intends  to  use  the  polynomial  Pzn^  1  (  ^  )  in  place 
of  the  function  *  ($  )  one  will  get  from  equation  (2) 

fl  - 

-  a  |(£* )  %  (• 5s)  + 

-/ 

But  equation  (10)  had  proved  the  integral  to  be  independent  of  the 
polynomial  %  (§ >■  hence  one  is  oblised  to  take  the  vaiaes  °f  ^ 
which  are  the  roots  of  the  Legendre  polynomial  Pn*^j  (  because  at 

these  points  this  function  vanishes. 


(12) 


Suppose  one  defines  another  polynomial 


- (?~*W 


(13) 


then  an  expansion  of  the  fraction 
fractions  yields 


+ 1  (§>  )  into  partial 


±\ 


21  gsiS) 

*  fvvn  C?* )  (H&) 

41 


r =  l/ 

2l  *  K-i  FnnGS, )($-*») 

Then  from  formulae  (ll)  and  12) 


(14) 


(15) 


(39) 


The  weights  W  ^  are  given  by 
+< 


— I 

and  comparison  with^  (15)  giv^e-s  the  values 


= 


F-n-H  C  %  )  dtp 


The  values  of  and  ^  for  Gaussian  quaduature  formulae  up  to  n 

(17) 


=  16 


have  been  tabulated  by  Lowan,  Davids  and  Levenson  and  were  the  ones 
used  during  the  present  investigation. 

The  method  of  energy  minimization  was  carried  out  in 


two  parts: 

(l)  Minimization  of  the  linear  parameters:- 

This  is  essentially  a  solution  of  the  matrix  equation  of 


the  type:- 

(  H  -  ES)  C  *  0  (17) 

Where  H  is  the  matrix  representation  of  the  Hamiltonian  operator  under 

the  chosen  set  of  basis  functions,  S  is  a  matrix  representation  of  1,  the  id- 

■> 

entity  operator,  under  the  same  set,  C  is  a  vector,  the  components  of  which 

are  the  linear  parameters  to  be  optimized,  and  E  is  the  eigenvalue 

— ^ 

corresponding  to  a  given  C  vector.  The  minimization  procedure  then 

— * 

essentially  consists  in  the  generation  of  a  C  vector  which  will  give  the 
lowest  E  value  which  might  be  said  to  correspond  to  the  ground  state  energy. 

An  examination  of  equation  (17)  shows,  that  it  can  hold  either 
if  the  vector  C  or  the  determinant  is  zero 

)  J±  —  fS  )  =  o 


(18) 


- 


(40) 


the  former  of  course  corresponds  to  an  uninteresting  trival  case  and  hence, 
one  has  to  look  for  vanishing  determinant.  It  is  obvious  then  that  (18)  is 
a  polynomial  h  :aving  several  roots  corresponding  to  each  there  being  ». 
separate  C  vector,  and  hence  to  generalize  even  further  one  may  replace 

— V  — V 

the  C  vector  by  a  matrix  C  compounded  from  several  vectors/  C  each 

corresponding  to  a  given  value  of  E 

(H  ES)  C_-  0  (19) 

where  E  also  takes  on  matrix  form,  being  diagonalized  with  the  diagonal 

elements,  equal  to  the  eigenvalues. 

The  method  used  for  the  solution  of  the  matrix  equation 

was  the  Jacobi  method;  but  unfortunately  this  cannot  handle  a  problem  where 

the  S  matrix  is  not  a  unit  matrix.  Let  us  then  examine  a  sligtt^different 

problem . 


yr\  -  f  (ta  s in)  > 


(20) 


Where  m  and  ^  are  new  matricies,  the  commutation  of  E  and  m  being 
permissible  since  the  former  is  a  diagonal  matrix.  The  matrix  m  is  So 

This  process  of  reduction  to  a  unit  matrix  is  known 


chosen  that 


- 


I 


4 


(41) 


V 


as  the  Schmidt  orthogonalization  process.  It  is  a  stepwise 
conversion  of  a  set  of  non-  orthonormalized  functions  into  an 
orthonormal  b&si®  set.  Suppose  one  has  two  non-orthonormal 
functions ,  and  X^  such  that 

Oi]*\>  41 

then  one  can  f  orm  by  taking  appropriate  linear  combination  of 
these;  two  other  functions  T^^and  7^,  in  the  .manner 


ts.  ax, 

(22) 

9^  ~  CX|4  ct 

(23) 

S 

<*/7  X/>  H 

(24) 

<*/  1  xi.>  ^  o 

(25) 

(26) 

This  then  gives  one  thjj-ee  unknown  linear  parameters,  a,  c,  and 
d  and  three  simultaneous  equations  (24),  (25)  and  (26)  from 
which  they  can  be  readily  evaluated.  If  one  writes  the  original 


S  matrix  from  the  non-orthonormalized  set  in  the  form 


■ 

' 

■  •  ■ 

' 

, 

•  ■  • 


(«) 


s  =  fa  r  SfcN  <*/!*■>  \ 

\Sa\  S «/  \  <&!*> 

one  has  essentially  by  a  Schmidt  orthogonalization,  convented 


into 


Xa->\  |R>I>  <(<**1 

<<:>,-* d>z  ja.X|  >  <^oti  +d j 


The  origin  of  this  matrix  beoormes  clear,  by  referring  to 
Dirac's  dual  vector  space  system  of  bras  and  kets  in  the 

fo  llojving  matrix  equation 


LA  O 


«  "j  {CXrfd»> 

d /\<W)  d  /  \<ti<i-WJii|cW>4p(,4dX2/^fW^ 


Where  the  matrix  on  the  left  has  been  multiplied  into  the  bra 
part,  while  the  one  on  the  right  into  the  ket  part  of  the  dot 
product.  These  two  left  and  right  matricides  are  the 
matricies  of  oqqation  (21). 


Now  multiplying  both  sides  of  (21)  by  on 


the  left 


r *  ^ 

iLTOA  —  >  E  >  >  A 


since  h  is  unitary  and 


(43) 


'  1*  =1. 

The  jA  matrix  is  so  chosen  that  the  matrix  product 

/-V/ 

/\  ]4  yy)  7v_ 

is  diagonalized,  then  the  resulting  diagonal  elements  will  be 
equal  element  by  element  to  the  matrix  E  and  hence  the 
eigenvalues  E  are  knowp,  onoethe  matrix  jd  H  id  has  been 
diagonalized. 

From  equation  (19)  one  has  again 

H  C  ESC 

It  follows  from  comparison  with  the  preceding  equation  that: 

C  :m  ^ 

Hence  the  vector  corresponding  to  the  lowest  diagonal  element 

— * 

in  E  wil  contain  the  best  set  of  C  values  for  a  miniminization 
of  energy. 

(2)  Minimization  of  the  exponential  parameters:-  : 

The  exponential  parameters  of  the  wave  functions 

were  minimized  quite  readily  by  a  method  of  Quaduatie  Minimiz- 

(18) 

ation  due  to  Ransil  ,  which  assumes  a  parabolic  energy  surface 


in  hyper  space:- 


■ 


(44) 


This  method  consists  of  effecting  a  MacLeurin 

expansion  of  the  energy  function  very  close  to  the  minimum. 

Considering  the  energy  to  be  a  function  of  two  exponential  variables 

only,  m,  and  m  one  can  construct  a  three  dimensional  plot, 

Z 

which  will  represent  the  energy  surface  for  a  given  functional 
form,  and  a  gi\*en  pair  of  numbers  m(  and  m^  will  determine 
the  minimum  of  this  surface.  Expressing  the  energy  as  a 
function  of  two  variables  in  the  form  E(m(,  m  )  a  MacL&urin 


expansion  about  the  zero  point  can  be  readily  obtained 


E(m,  »2>—EB+t(f^)#*n+£  i-X 


Then  applying  the  minimization  criterion  of 


Since  E  is  a  constant  and  i  and  j  are  dummy  indices  which  can 
be  readily  inter  changed,  one  can  simplify  the  above  equation  to 

«e-  qsy  «>«>+ 


(27) 


■  ' 


(43) 


This  equation  now  implies  that  either  all  are  zero  and 

no  incrementation  of  the  parameters  have  been  effected  or 
that  the  coefficients  of  (27  )  are  vanishing,  the  former  of  course 
corresponds  to  a  trivial  case,  and  the  problem  reduces  to  one  of 
the  type 


3*»t, 


a«j 


i-l,  ^ 


\ 

This  can  then  be  written  in  matrix  form: 


The  problem  is  then,  one  of  calcuhting  the  "best"  set  of 
incrementsJ2*Wfl(  and  /SnQ^as  shown  in  the  vector,  on  the  right 
hand  side  of  (28  ). 

For  a  given  set  of  mand  ^lvalues,  which  should 
be  very  close  to  the  minimum  one  proceeds  to  evaluate  the 
elements  of  the  matrix  in  the  above  equation  (24)  with  the  help 
of  Newton's  interpolation  formulae  for  one  and  two  variables; 
and  are  of  the  following  forms:- 

„  ££Cmr •  -  j'l 

v 


. 


The  inverse  of  the  matrix  is  evaluated  and  multiplied  into  the 
left  hand  vector,  resulting  there  by  in  an  improved  set  of 
parameters,  which  are  added  again  to  the  original  m  parameters 
and  an  optimized  energy  calculated. 

Programming  Procedures:- 


For  the  present  project  the  integrals  and  the 


various  matrix  manipulations  involved  were  programmed  for  the 
IBM  1620  made  available  by  the  University  of  Alberta  Computing 
Center.  A  general  flow  chart  showing  the  various  steps  in 
calculation  is  shown.  Those  routines  marked  with  an  asteriske, 
as  well  as  the  floating  point  subroutines  used  throughout  were 
programmed  by  F.  W.  Birss. 


Since  (n  +  l)  (  n-(-2)/2  energy  values,  were  required 


to  effect  the  <£uadrati'  minimization  procedure,  the  latter  was 
given  full  control  of  the  program. 


The  purpose  of  most  of  the  subroutines  is  given 


in  the  discussion  above.  The  eigenvalue  search  subroutine, 
was  designed  to  search  for  the  lowest  ei^ihvalues,  and  correspond¬ 
ing  eigenvectors,  and  to  set  these  in  preassigned  positions. 


! 


(47) 


Fig.  1:  Generalized  Flow  Chart  for  Program  used 


CHAPTER  111 


THE  FIRST  TRIAL  FUNCTION 


The  Trial  Function:- 


(48) 


In  order  to  compare  the  efficiency  of  the  Gaussian  approximation 


in  this  simplest  molecular  problem  with  the  normal  approach  using  Slater 
type  orbitals,  a  trial  function  comparable  to  that  of  Finkelstein  and  Horowitz, 
was  chosen.  This  had  the  form 


The  linear  parameters  (Ci)  and  the  exponential  parameters  (mi)  were 


obtained  by  the  method  given  in  Chapter  11.  The  results  obtained  are  given 
in  Tables  1  to  1 1 1 .  In  Table  1  V  are  given  the  results  of  a  calculation  of 
the  value  of  the  function  along  the  internuclear  axis  for  the  two  and  three 


Gaussian  combinations  at  Rc-Z.  In  figure  2  the  results  for  the  three  Gaussian 


"^ombinatiorB are  plotted  with  those  obtained  from  the  exact  and  from  the 


Finkelstein  -  Horowitz  treatment;  in  fig.  3  the  results  of  the  two  and  three 
Gaussian  Combinations  are  plotted  for  comparison.  . 

Discussion:- 


combination  of  Is  orbitals  centered  on  eachnucleus.  The  results  however 
were  not  very  accurate  as  no  allowance  had  been  made  for  an  orbital  exponent. 
The  lack  of  this  factor  results  in  a  decrease  in  kinetic  energy  of  the  molecule, 
instead  of  an  increase  as  might  be  expected  for  stable  bond  formation.  The 
'Quantum  Mechanical  Virial  theor/em  is  violated.  At  an  interP  pc  lear  separation, 
close  to  the  equilibrium  value,  the  influence  of  a  second  nucleus  is  likely 
to  cause  an  increase  in  the  amount  of  positive  charge  an  electron  would 


experience  as  a  result  of  being  in  a  molecular  orbital.  The  introduction  of 


y 


V  i 


(49) 


an  effective  nuclear  charge  at  both  ends  of  the  molecule  woulcj  result  in  an  in¬ 
creased  attraction  for  the  electron, intothe  region  between  the  nucle  ii,  and 
thereby  contribute  to  stable  bond  formation.  Minimization  of  the  energy  with 
respect  to  these  parameters  does  in  fact  render  the  wave  function  consistent 
with  the  virial  the  ore  in  in  keeping  with  an  important  theoretn,  proved  for 
the  first  time  by  Fock^^  . 

As  the  initial  trial  value  of  the  exponential  parameters 

in  the  Gaussian  function,  the  value  of  the  parameters  obtained  by  Longstaff 

(23)  A 

and  Singer  for  the  hydrogen  atom  were  scaled  to  the  hT?  problem  by 

multiplying  by  the  square  of  the  orbital  exponent  obtained  by  Finkelstein  and 
Horowitz  who  had  for  the  first  time  realized  the  necessity  for  such  a  para¬ 
meter.  The  result  for  the  electronic  energy  at  internuclear  distance  of 
two  atomic  units  was  found  to  be  -T  0851  a.  u,  the  value  of  the  parameters 
being  n^«0‘7179  =0*3393  rr5*O“0358  for  a  three  Gaussian  combination. 

From  an  inspection  of  Table  1  it  is  obvious  that  a  considerable  variation  of 
the  exponential  parameters  was  effected  to  obtain  the  minimum  value  of  - 
1-0882  a.  u. 

The  energies  obtained  for  various  linear  combinations  of 
this  type  of  function,  do  indeed  compare  very  favorably  with  the  results  of 
Finkelstein  and  Hoitwitz.  However  the  value  of  the  function  as  plotted  against 
distances  along  the  internuclear  axis  do  not  compare  at  all  favorably  with 


the  exact  or  the  Finkelstein  -  Horowitz  curve. 
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The  character  of  the  values  of  the  wave  function  along 
the  internuclear  axis  indicate  very  clearly  that  the  Gaussian  approximation, 
although  leading  to  a  reasonable  value  for  the  energy,  cannot  reproduce 
in  detail  the  electronic  distribution  curve.  The  shape  of  the  function  near 
the  nucleii  still  exhibits  a  considerable  amount  of  the  Gaussian  character, 
in  that  the  tangent  at  the  maxima  appear  to  be  zero. 

Examination  of  Fig  3  shows  some  of  the  nature  of  the 
improvement  to  be  expected  as  one  increases  the  number  of  basis  functions. 
Although  the  electronic  distribution  in  the  center  of  the  bond  is  not  improved, 
the  character  of  the  function  near  the  nucleii  and  outside  the  nucleii,  show 
marked  improvements.  The  two  Gaussian  functiorC  exhibits  a  rather  flat 
maximum  which  is  distinctly  displaced  into  the  center  of  the  bond  (This 

(24") 

shift  of  the  maxima  was  in  fact  introduc  ed  by  Gurnee  and  Magee'  -for  the 
hydrogen  molecule,  who  found  a  displacement  of  0°11  a.  u  toward  the  center  of 
the  molecule).  On  the  other  hand  the  three  Gaussian  combination  displays 
a  rather  sharper  maxima  which  is  more  towards  the  character  of  the  exact 
solution.  Also  it  falls  off  more  sharply  than  the  two  Gaussian  in  the  region  out  = 
side  the  nucleii.  One  therefore  concludes  that  the  improvement  in  the 
energy  may  be  attributed  to  the  simultaneous  improvement  of  the  function 
in  the  immediate  region  of  the  nucleii. 

With  a  number  of  rather  drastic  assumptions  one  can  attempt 
an  extension  of  the  present  results  to  higher  homonuclear  diatomic  molecules. 
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If  one  assumes  that  the  character  of  the  ls£  orbital  is  not  much  affected 
by  double  occupancy  or  by  the  presence  of  electrons  in  other  orbitals,  then 
one  may  say  that  the  energy  of  the  ls4  orbital  will  vary  roughly  as  the  square 
of  the  charge  on  the  nucleii.  Hence  the  percentage  error  introduced  by 
the  Gaussian  approximation  for  the  ls<£  orbital  will  remain  roughly  constant. 
The  absolute  value  of  this  discrepancy  will,  however,  increase.  Although 
this  may  not  have  a  drastic  relative  effect  upon  the  total  electronic  energy, 
the  error  so  introduced  into  bond  dissociation  energy  predictions  (which 
remains  very  roughly  constant  from  molecule  to  molecule)  might  be  fairly 
high.  From  this  qualitative  argument  one  might  expect  that  the  Gaussian 
approximation  for  lsdorbitals  in  higher  molecules  may  lead  to  rather  large 
errors  in  the  prediction  of  bond  dissociation  energies. 

In  the  application  of  Gaussians  to  multi  electronic  problems  , 
this  unfortunate  feature  may  be  com  pensated  by  the  use  of  a  larger  number 
of  Gaussian  functionsper  Slater  orbital.  However  there  must  be  made  an. 
attempt  to  limit  the  number  of  basis  functions  so  chosen  in  order  to  remain 
within  the  practical  limits  imposed  by  the  demands  of  computation. 

Against  this  unfavorable  feature  one  may  set  the  advantage  of 
Gaussian  functions  in  their  application  to  explicit  electronic  correlation. 

.r 

j 

functions.  Further  investigation  into  multi  electronic  problems  are  needed 
to  give  some  indication  whether  this  feature  may  sufficiently  compensate 
the  use  of  Gaussian  functions  as  compared  to  Slater  orbitals. 


f 


'toq  . 


. 


Table  1 
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! 


i 

Results 

from  the  first  trial  function:  R  : 

2  3  4 

=  2  a.  u.  * 

5 

6 

E 

e 

(a.  a.  ) 

-1. 0803 

-1. 0882 

-1. 0888 

-1. 0906 

-1.0909 

E 

m 

(a.  a.  ) 

-0. 5803 

-0. 5882 

-0. 5888 

-0. 5906 

-0. 5909 

ml 

2.0056 

0. 2803 

0.  1725 

0.  1947 

0. 1392 

m2 

0. 3343 

1. 1687 

0. 6582 

0 . 4666 

0. 3891 

m  ^ 

- 

8. 1750 

0.  3014 

1.4658 

0. 9469 

m4 

- 

- 

17. 1689 

5. 5351 

2. 3044 

rr>5 

- 

- 

- 

34. 4770 

8. 8075 

m6 

- 

- 

- 

- 

56. 9970 

c  i 

0. 1972 

0. 7429 

0. 4515 

0. 3575 

0. 2209 

C  2 

0. 8653 

0.2974 

0. 5505 

0. 5207 

0. 6034 

C  3 

- 

0. 0400 

0. 0848 

0. 1598 

0. 1632 

C4 

- 

- 

0.0157 

0. 0435 

0. 0887 

C5 

- 

- 

- 

0. 0068 

0. 025^ 

C6 

- 

- 

- 

0. 0036 

* 

In  this 

and  the  following  tables 

R  indicates  internuclear 

separation, 

Ee  the 

calculated  electronic  energy  and 

E  rn  the  calculated  molecula 

energy,  i  indicates 

the  numbe 

r  of  basis 

functions  used. 

. 
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Table  1 1 


Results 

from  the  fir 

st  trial  function:  R  = 

1 . 8  and  2.2  a.  u. 

R  (a.  u.  ) 

1.8 

1.8 

2.2 

2.2 

i' 

3 

6 

3 

6 

Ee(a-  u.  ) 

-1.1415 

-1.  1444 

-1.  0411 

-1. 0434 

Em(a.  u.  ) 

-0. 5860 

-0. 5888 

-0. 5865 

-0.  5889 

ml 

0*2857 

0. 1405 

0.2731 

0.  1392 

m  2 

1. 1687 

0. 3944 

1.  1687 

0. 3893 

m  3 

8. 3420 

0. 9540 

7. 9860 

0. 9438 

m  4 

- 

2. 3298 

- 

2. 3052 

m  5 

- 

8. 8634 

- 

8. 7666 

m  6 

- 

5. 6880 

- 

56. 6568 

C  1 

0. 7088 

0. 1887 

9. 7680 

0. 2425 

C  2 

0.  3312 

0. 5843 

0. 2735 

0.  6041 

C  3 

0. 0411 

0. 2027 

0. 0395 

0.  1306 

C  4 

- 

0. 0863 

- 

0.  0914 

C  5 

- 

0. 0275 

- 

0. 0241 

c  6 

- 

0. 0038 

- 

0. 0036 

'  - 
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Table  1 1  1 


Results  from  the  first  trial  function:  three  Gaussian  combination 


R  (a.  u.  ) 

10  ~29 

1.9 

2.  1 

o 

o 

fM 

E  e  (a.  u.  ) 

-1. 9702 

-1. 140 

-1. 0639 

-0. 5443 

E  m  (a.  u.  ) 

-1. 9702 

-0. 5877 

-0. 5877 

-0. 4943 

m  l 

0. 2970 

0. 2832 

0. 2769 

0. 1883 

m2 

1.  1687 

1. 1685 

1. 1688 

1. 1681 

m  3 

0.  9413 

8. 2610 

8. 0870 

9.  9291 

C  1 

0.2753 

0. 7270 

0. 7565 

0. 7930 

C  2 

0.  6626 

0. 3130 

0. 2843 

0. 2946 

C  3 

0.  1991 

0.  0405 

0. 0397 

0. 0291 
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Table  IV 


Value  of  wave  function  along  intern  uclear  axis 


r  A  (a.  u.  ) 

r  B  (a*  u-  ) 

i=2 

i-3 

1.0 

3.  0 

0.  142 

0.  050 

0.  3 

2.3 

0.  317 

0.  330 

0.2 

2.2 

0.  338 

0.  368 

0.  1 

2.  1 

0.  354 

0.  398 

0.0 

2.0 

0.  363 

0.  412 

0.  1 

1.9 

0.  365 

0.  408 

0.2 

1.8 

0.  360 

0.  387 

0.3 

1.  7 

0.  350 

0.  355 

0.  8 

1.2 

0.  278 

0.279 

0.  9 

1.  1 

0.  270 

0.  275 

1.0 

1.0 

0. 26& 

0.  274 

CHAPTER  IV 


THE  SECOND  TRIAL  FUNCTION 
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The  Trial  Function:- 

One  of  the  classic  ways  of  improving  the  electronic  distribution 

function  in  the  hydrogen  molecular  ion  problem  has  b^en  to  introduce  into 

the  trial  function  a  contribution  whose  main  weight  lies  along  the  inter- 

(9) 

nuclear  axis.  Walsh  and  Matsen  ,  from  a  one  dimensional  treatment 
showed  the  importance  of  considering  the  polarizability  of  the  electron 
cloud,  in  fact  these  workers  used  a  hundred  percent  constraint  to  the  inter- 
nuclear  axis  and  obtained  a  marked  improvement  in  the  binding  energy. 
Dickenson  introduced  the  Slater  orbital  into  the  trial  function  and 
found  an  improvement  in  the  energy,  along  with  a  satisfactory  improvement 
in  the  axial  distribution  curve. 

These  arguments  prompted  the  choice  of  the  second  function 

y._  t-CL  *&-**■**>-*$  ?***> 

Which  involves  the  introduction  of  a  weighting  parameter  q  along  the  Z-axis, 
chosen  to  be  the  internuclear  axis. 


Piscussion:- 


It  is  recognized  that  this  form  is  not  the  most  satisfactory, 

although  it  is  the  only  one  which  maintains  the  simple  Gaussian  character 

< 

throughout  the  trial  function  along  with  some  aspects  of  polarization.  Its 
deficiency  arises  from  the  impartial  character  of  Z  in  that  the  electron 
distribution  along  the  internuclear  axis  is  effected  equally  both  between 
and  outside  the  nucleii.  This  does  result  in  a  relative  increase  in  the 
electron  charge  density,  in  the  region  of  the  nucleii;  but  contrary  to  the 
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function  of  Dickenson  (and  that  described  in  Chapter  V)  this  does  not  arise 
from  a  depletion  of  electron  charge  density  in  the  region  beyond  the  nucleii. 
Rather  any  building  up  of  charge  density  in  the  internuclear  region  must 
be  effected  at  the  expense  of  the  relative  charge  density  in  the  region  of  the 
nucleii. 

That  this  is  not  an  entirely  favorable  feature  is  pointed 
out  by  a  comparison  of  the  results  of  this  and  previous  Chapter  (using-  the 
three  Gaussian  function  as  an  example).  There  is  a  relatively  minor  effect 
upon  the  energy  -1*0882  a.  u.  vs  -1-0894  a.  u.  for  the  first  and  second  trial 
functions  respectively.  Further  one  detects  a  tendency  for  resistance 
on  the  part  on  the  function  for  any  great  deviation  from  spherical  character. 

In  view  of  the  minor  improvement  in  the  energy,  and  of  the 
amount  of  computation  involved  in  the  necessary  numerical  integration 
the  conclusion  is  that  such  an  attempt  to  weight  the  distribution  along  the 


internuclear  axis  is  not  recommended. 
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TABLE  V 


Results  from  the  second  trial  function:  R  =  2  a.  u. 


i 

2 

3 

4 

E  e  (a-  u-  ) 

-1. 0785 

-1. 0894 

-1.0909 

(a.  u.  ) 

-0. 5785 

-0. 5894 

-0. 5909 

ml 

2. 0056 

0. 2630 

0. 1823 

m  2 

0. 3037 

1.  2177 

0. 6582 

m 

- 

9.5299 

2, 8858 

m  4 

- 

- 

19. 1696 

^  1 

2. 0056 

0.2997 

0. 2278 

<1  2 

0. 3037 

1.  0776 

0. 6582 

q  3 

- 

8. 5037 

2. 8858 

q  4 

- 

- 

9. 1695 

C  1 

0. 2083 

0. 7846 

0. 4769 

C  2 

0. 8082 

0. 2777 

0.  5141 

C  3 

- 

0. 0435 

0,09209 

C  4 

0. 0147 
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TABLE  VI 


Results  from  the  second  trial  function:  three  Gaussian  combintation 


R  (a.  u.  ) 

1.8 

2.  1 

2.2 

20 

E  (a.  u.  ) 
e 

-1. 1424 

-1. 0631 

-1. 0414 

-0. 5362 

E  (a.  u.  ) 
m 

-0. 5868 

-0. 5869 

-0.  5869 

-0. 4862 

m 

1 

0. 2736 

0. 5453 

0. 2530 

0.  2219 

mz 

1.  2376 

2. 5561 

1.2562 

1. 3208 

m 

3 

9. 1642 

23. 8021 

7. 7280 

8.  5225 

qi 

0.2997 

0.2997 

0.2997 

0.2997 

q2 

1.  1481 

1. 8496 

1. 1479 

1.  1907 

q  3 

8. 5623 

5. 5674 

8. 7238 

8. 6451 

C 

1 

0. 1993 

0. 7990 

0. 8112 

0. 8376 

c  z 

0. 5457 

0. 2655 

0. 2554 

0. 2751 

C  3 

0.2202 

0. 0431 

0.  042  9 

0. 0316 

CHAPTER  V 


The  Third  Trial  Function 
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The  Trial  Function:- 

We  had  seen  in  Chapter  111,  how  the  idea  of  Finkelstein 
and  Horowitz  of  introducing  an  effective  nuclear  charge  in  the  case  of  the 
Slater  orbitals  had  served  to  effect  &  marked  improvement  in  the  energy, 
one  of  the  reasons  being  the  presence  of  a  larger  nuclear  charge, 

resulting  in  a  greater  attraction  for  the  electron,  causing  thereby  an  increased 
electronic  concentration  in  the  internuclear  region.  Following  these  lines 
of  argument,  and  the  work  of  Walsh  and  Matsen,  we  tried  the  second  type 
of  function  with  the  idea  of  weighting  along  the  Z  axis,  chosen  to  be  the 
co-ordinate  along  which  the  internuclear  axis  lies.  However  this  resulted 
in  an  increase  of  the  electronic  density  in  the  middle  of  the  two  nucleii,  at 
the  expense  of  that  in  nuclear  regions. 


Dickenson  improved  upon  the  Finkelstein  -  Horowitz  result 


by  the  introduction  of  a  2p  orbital.  This  orbital  when  centered  upon  the 
nucleii  would  make  a  negative  contribution  towards  the  electronic  cloud  in 
the  outer  regions,  and  a  positive  contribution  in  between  the  two  nucleii. 
The  introduction  of  this  p  orbital  did  in  fact  markedly  improve  the  energy. 


which  in  effect  is  a  linear  conbination  of  two  s  type  functions  and  one  p  type 


function,  centered  at  each  nucleus. 


« 
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Discussion:- 

At  an  internuclear  separation  of  two  atomic  units  a  molecular 
energy  of  -0.5878  a.  u.  was  obtained,  as  against  values  of  -0.5803  a.  u.  and 
-0. 58^2  a.  u.  from  linear  conbination  of  two  and  three  s  type  Gaussian 
functions  respectively. 

This  alone  suggest  that  with  the  restricted  practical  basis 
set  used  the  problem  is  not  necessarily  '  that  of  increasing  the  electron 
distribution  in  the  internuclear  region  but  an  improvement  over  the  entire 
molecule,  including  the  region  at  the  nucleii. 

An  examination  of  the  results  given  in  Tables  Vlll  and  figures  3 
and  4  shows  that  there  is  minimal  improvement  in  the  electronic  distribution 
over  the  whole  molecule  as  one  increases  the  basis  set  from  two  s  type 
Gaussian  functions  to  the  combination  of  the  latter  and  one  p  type  function. 

It  is  indeed  mor«  favorable  to  use  s  type  function  as  a  third  member  of 
the  basis  set. 

This  result  does  not  parallel  that  of  Dickenson  who  found  an 
improvement  of  0.0175  a.  u.  in  the  energy.  The  prime  difference  between 
the  basis  set  used  by  Dickenson  and  that  used  here  is  the  Gaussian  form  of 
the  latter.  In  both  cases,  however,  the  approximate  contribution  of  the 
p  type  wave  function  is  under  5%.  This  would  seem  to  indicate  that  the 
expected  effect  of  the  p  type  function  is  being  offset  by  the  fact  that  the 
small  number  of  s  type  Gaussian  functions  used  do  not  accurately  describe 
the  distribution  near  the  nucleii.  A  Is  Slater  function  properly  describes 


(62) 


certain  of  the  characters  of  the  exact  function  in  this  region,  in  particular 
the  discontinuities-  in  the  tangent  at  the  nucleii.  The  results  given  in  Chapter 
111  show  that  thes-e  characteristics  can  be  approached  by  increasing  the  number 
of  s  type  functions- in  the  basis  set.  The  conclusion  to  be  drawn  from  the 
above  results  are  that  such  an  increase  is  necessary  before  one  can  expect 
an  independant  improvement  by  the  introduction  of  a  p  type  function.  An 
enlarged  basis  set  would  however,  imply  evaluation  of  a  very  large  number 
of  matrix  elements-  and  of  individual  energies  for  an  optimization.  The 
amount  of  computation  required  for  both  of  these  portions  depends  approximately 
on  the  second  power  of  the  number  of  functions  in  the  basis  set. 

From  a  practical  point  of  view  one  may  say  that  the  application 
of  Gaussian  functions  to  molecular  orbital  calculations  is  limited  by  the 
availibility  of  high  speed  computers  of  large  t  caipacity.  This  conclusion 
however  does  not  exclude  the  possible  use  of  a  small  basis  set  of  Gaussian 
functions  in  problems  involving  the  explicit  use  of  inter  electronic  co¬ 
ordinates  where  Slater  functions  introduce  rather  formidable  integrals. 
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Table  VI 1 


Results  from  third  trial  function:  Linear  combination  of  three  Gaussian 


functions 

R  a.  u. 

1.8 

2.0 

2.2 

2.4 

E 

e 

-1.  1401 

-1. 0878 

-1. 0417 

-1. 0023 

E 

rn 

-0. 5S45 

-0. 5878 

-0. 5871 

-0. 5857 

ml 

2. 3342 

2.1342 

2.3046 

2. 3915 

m2 

0. 3831 

0. 3471 

0. 3471 

0. 3486 

3 

0.08Q02 

0. 1001 

0. 0790 

0. 0791 

C1 

0.  12.34 

0. 1296 

0. 1165 

0. 1097 

C2 

0.  4643 

0.4653 

0. 4859 

0. 5001 

c  , 

0. 1257 

0. 1086 

0. 1201 

0.  1225 
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Table  VI 11 


Value  of  function  along  internuclear  axis. 


r  A 

r  B 

s  contribution 

p  contribution 

Total  (p  functi 

1.0 

2.0 

0.  1313 

0. 0027 

0. 1340 

0.  3 

1.  3 

0. 3027 

0. 0092 

0. 3119 

0.2 

1.2 

0. 3244 

0.0101 

0. 3345 

0.  1 

2.  1 

0. 3402 

0. 0109 

0. 3511 

0.0 

.1.0 

0. 3492 

0.  0117 

0. 3609 

0.  1 

0.9 

0. 3507 

0.0124 

0. 3631 

0.2 

0.8 

0. 3454 

0.0131 

0. 3585 

0.  3 

0.  7 

0. 3341 

0.0137 

0. 3478 

0.8 

0.2 

0. 2594 

0. 0156 

0. 2750 

0.9 

0.1 

0.  2.538 

0.0157 

0. 2695 

1.0 

0.0 

0. 2535 

0. 0158 

0. 2693 
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Table  IX 

Comparison  of  energy  at  internuclear  distance  of  2  a.  u 


Burrau 

-0.  5987 

Hyelleraas 

-0. 5993 

Jaffe 

-0. 5993 

Guillimin  and  Zener 

-0. 599^ 

Finkelstein  and  Horowitz 

-0. 5799 

Dickenson 

-0.5976 

Present  work:- 


First  trial  function 

i  -  2 

-0. 5803 

i  =  3 

-0. 5882 

i  -  4 

-0. 5888 

i  =  5 

-0. 5906 

i  -  6 

-0. 5904 

Second  trial  function 


i  =  2 

-0. 5785 

i  =  3 

-0. 5899 

i  --  4 

-0. 5909 

Third  trial  function 


i  =  3 

-0. 5878 

Chapter  VI 


Spectroscopic  Constants 
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In  Tabled  are  collected  the  data  used  in  calculating  the 

various  spectroscopic  constants.  The  constants  are  displayed  in  TableWj. 

The  formulae  used  in  this  calculation  were  taken  from  a  publication  by  Frag  a 

(25) 

and  Ransil 

The  majority  of  the  results  are  surprisingly  good  considering 
the  comments  which  have  been  made  in  the  preceding  Chapter-.  Fraga  and 
^Ransil  from  a  study  of  a  wide  variety  of  diatomic  molecules  conclude  The 
rather  regular  behavior  suggests  the  possibility  ofu,y  and^/,as  criteria  for 
the  ^goodness  of  a  wave  functiorf^  .  although  they  detect  no  correlation 
between  the  accuracy  of  these  quantities  and  that  of  the  calculated  total 
energies  or  dissociation  energies. 

The  present  results,  coupled  with  the  information  obtained 
on  the  actual  values  of  the  wave  function,  cast  severe  doubt  upon  any  safe 
prediction  of  the  goodness  of  the  wave  function  based  on  derived  spectro¬ 
scopic  constants.  In  this  regard  it  may  be  reported  that  a  preliminary  invest¬ 
igation  using  an  assuredly  absurd  wave  function  gave  rise  to  spectroscopic 
constants  which  were  as  good  as  or  better  than  reported  by  Fraga  and  Ransil. 


r 
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Table  X 

Data  used  in  calculation  of  spectroscopic  constants 


First  function 


R  =  1.8 

E  m"  -0 . 58597 

1.9 

»  -0..5877.L 

2.0 

-  -0.58821 

2.  1 

-  -0.58775 

Third  function 

R  -  1.  8  E  m=  -0.  58454 
.2.0  -  -0.58776 

2.2  =  -0.58715 

2.4  •  -0.58570 


Table  XI 

Spectroscopic  Constants 

_  5 

^cmH^cmjBeCcm'1)  ^cm"1)  Re(A°)  Em(a-U-)  10”  ke(dynes/cm) 


First  function 

2249 

197 

30.  0 

2.  1 

1. 057 

-0.  5.882 

1.  501 

Third  function 

2024 

554 

28.  5 

5.4 

1. 083 

-0. 5879 

1.216 

(25) 

Fraga  and  Ransil 

2275 

227 

28.  0 

2.2 

1. 093 

-0. 5791 

1. 536 

Experimental 

2297 

62 

29.8 

1.4 

1.  06 

-0. 6026 

1.  566 
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Distance  from  the  center  of  bond  (a. 
Fig.  2:  Axial  plot  of  wave  function 


Value  of  wave  function 
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Axial  plot  of  the  wave  function 


Value  of  wave  function 
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Distance  from  the  center  of  bond  (a. 
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